Power laws and stretched exponentials in a noisy finite-time-singularity model 
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I. INTRODUCTION 

There is a continuing current interest in the influence 
of noise on the behavior of nonlinear dynamical systems 
jl], ||] . Here the issues for example associated with the in- 
terference of stochastic noise with the deterministic chaos 
of maps H or extended systems, as for example the noise- 
driven Kuramoto-Shivashinski equation [Q, are of funda- 
mental interest. 

In a particular class of systems the nonlinear charac- 
ter gives rise to finite-time-singularities, that is solutions 
which cease to be valid beyond a particular finite time 
span. One encounters finite-time-singularities in stellar 
structure, turbulent flow, and bacterial growth [p|, p\ |Tf a s 
well as in in Euler flows and in free-surface- flows ]8|, S, . 
Finitc-time-singularities are also encountered in model- 
ing in econophysics, geophysics, and material physics 

In the context of hydrodynamical flow on a nanoscale 
[ p7t , where microscopic degrees of freedom come into 
play, it is a relevant issue how noise influences the hydro- 
dynamical behavior near a finite time singularity. Leav- 
ing aside the issue of the detailed reduction of the hy- 
drodynamical equations to a nanoscale and the influence 
of noise on this scale to further study, we assume in the 
present context that a single variable or "reaction co- 
ordinate" effectively captures the interplay between the 
singularity and the noise. 

We thus propose to consider a simple generic model 
system with one degree of freedom governed by a nonlin- 
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ear Langevin equation driven by Gaussian white noise, 

dx A 
~dt = ~2\x\ 1 +i> 



+ 77, (rpi)(t) = A8(t) . (1.1) 



The model is characterized by the coupling parameter A, 
determining the amplitude of the singular term, the index 
/i > 0, characterizing the nature of the singularity, and 
the noise parameter A determining the strength of the 
noise correlations. Specifically, in the case of a thermal 
environment at temperature T the noise strength A cx T. 

In the absence of noise this model exhibits a finite-time 
singularity at a time to, where the variable x vanishes 
with a power law behavior determined by fi. When noise 
is added the finite-time-singularity event at to becomes 
a statistical event and is conveniently characterized by a 
first-passage-time distribution W(t) For zero noise 

we thus have W(i) — S(t — to), restating the presence of 
the finite-time-singularity. In the presence of noise W(t) 
develops a peak about t — to, vanishes at short times, 
and acquires a long time tai l. 

The model in Eq. (1.1) has also been studied in 
the context of persistence distributions related to the 
nonequilibrium critical dynamics of the two-dimensional 
XY model jl9| and in the context of non-Gaussian 
Markov processes [^0|. Finally, regularized for small x, 
the model enters in connection with an analysis of long- 
range correlated stationary processes . 

It follows from our analysis below that for /i = 0, the 
logarithmic case, the distribution at long times is given 
by the power law behavior 



W(t) ~ r 



2A 



(1.2) 



For vanishing nonlinearity, i.e., A = 0, the finitc-timc- 
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singularity is absent and the Langevin equation (1.1) de- 
scribes a simple random walk of the reaction coordinate, 



2 



yielding the well-known exponent a = 3/2 [[l8], p3| . 
In the nonlinear case with a finite-time-singularity the 
exponent attains a nonuniversal correction depending on 
the ratio of the nonlinear strength to the strength of the 
noise; for a thermal environment the correction is pro- 
portional to 1/T. 

In the generic case for \i > considered here we find 
that the fall off is slower and that the correction to the 
random walk result is given by a stretched exponential 



W(t) ~ t~ 3/2 exp[A(t 



-m/(2+m) 



1] 



(1.3) 



where A — > A/A/i for fi 



0; 



i n th e limit fj, — > this 
expression reduces to expression (1.2). The above results 
for pk = and /i > have also been obtained in the 
context of the critical dynamics of the XY model Jn|] . 

The paper is organized in the following manner. In 
Section || we introduce the finite-time-singularity model 
in the noiseless case and discuss its properties. In Section 
[II we consider the noisy case and discuss the ensuing 
Langevin equation and associated Fokkcr Plank equa- 
tion. We discuss the relationship to an absorbing state 
problem and introduce the first-passage-time or absorb- 
ing state distribution. In Section ^ we review the weak 
noise WKB phase space approach to the Fokker-Planck 
equation, apply it to the finite-time-singularity problem, 
and discuss the associated dynamical phase space prob- 
lem. In Section |y| we apply the WKB phase space ap- 
proach and evaluate the weak noise absorbing state distri- 
bution at long times. We derive the random walk result 
in the linear case for A = 0, the power law tail for fi = 0, 
and the stretched exponential behavior for /i > 0. In Sec- 
tion VI we derive an exact solution of the Fokker-Planck 
equation in the case /x = in terms of a Bessel func- 
tion and present an expression for the absorbing state 
distribution. In Section |VI]| we present a summary and 
a conclusion. Details of the phase space method is dis- 
cussed in Appendix |a|; aspects of the exact solution in 
Appendix |B. 



II. FINITE-TIME-SINGULARITY MODEL 



Let us first consider the noiseless case for A = 0. It is 



instructive to express the equation of motion (1.1) in the 
form 



dx 
lit 



IdF 
2~dx 



where the potential or free energy has the form, 
F(x) — A log \x\ for /x = 0, 
F(x)=—\x\- fi for ^>0. 



(2.1) 

(2.2) 
(2.3) 



The free energy has a logarithmic sink for /x = and a 
power law sink for general fi > 0. In both cases F drives 
x to the absorbing state x = 0. Solving Eq. (2.1) in the 
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FIG. 1: In a) we show the time evolution of the single degree 
of freedom x. x reaches the absorbing state x — at a finite 
time to- In b) we depict the free energy F(x) driving the 
equation. The absorbing state x — corresponds to the sink 
in F(x). 



logarithmic case we obtain for positive x the solution 

x = V\Vt^~~t . (2.4) 

This solution displays a finite-time-singularity at to = 
Xq/X, where xq is the initial value at time t = 0, with 
x approaching the absorbing state with exponent 1/2. 
In other words, the attraction to the sink in the free 
ener gy o ccurs in a finite time span; for times beyond to 



Eq. ( 2.1 ) does not possess a real solution. This is the way 
we define a finite-time-singularity in the present context. 
For general /x > we obtain the generalization of the 
solution (EDI), 



x = [A(2 + A i)/2] 1 /( 2+ ^[t - i]V(2+M) 



(2.5) 



which approach the absorbing state with exponent 1/(2 + 
fi) at time t Q = 2xq +m /A(2 + /i). In Fig. [I] we have for 
fi = depicted the solution x and the free energy F(x) 
driving x to zero j52j . 



III. LANGEVIN AND FOKKER-PLANCK 
EQUATIONS 

In the presence of noise the finite-time-singul arit y 
problem is characterized by the Langevin equation ( |l . l| ) . 
Here the noise drives the variable x into a fluctuating 
state in competition with the free energy which tends to 
drive x to the absorbing state x = 0. In the case fi > 
treated here with an absorbing state the free energy has 
a sink and there is no stationary distribution; the prob- 
ability leaks out at x — Q . 

From another point of view, introducing the variable 



y 



the Langevin equation (|DJ) takes the form 



1 dy 
2 + fidt 



(3.1) 



In the noiseless case the variable y decreases linearly to 
zero at time t = t . In the presence of noise y fluctuates. 
We note, however, that the noise is manifestly quenched 
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at y = 0, yielding the absorb ing state. Absorbing state 
models of the type in Eq. (3.1) for extended systems 
have been studied extensively in the context of directed 
percolation, catalysis, and Reggeon field theory [E4l \28, 



In order to analyze the stochastic aspects of finite-time- 
singularities in the presence of noise we need the time de- 
pendent probability distribution P(x, t) and the derived 
first-passage-time or absorbing state probability distribu- 
tion W(t). The distribution P(x,t) is defined according 
toil El 



P(y,t) = (6(y - x(t))) 



(3.2) 



where x is a stochastic solution of Eq. (1.1) and (■ • •) 
indicates an average over the noise r\ driving x. In the 
absence of noise P(y,t) — S(y — x(t)), where x is the 
deter min istic solution of Eq. ( |2.l|) given by Eqs. (2.4) 
and (2.5) and depicted in Fig 111 At time t = the 
variable x evolves from the initial condition :ro, implying 
the boundary condition 



P(x,0) = S(x-x ) 



(3.3) 



At short times x is close to Xq and the singular term is not 
yet operational. In this regime we then obtain ordinary 
random walk with the Gaussian distribution 



P(x,t) = (2 7 rAi)- 1 / 2 



exp 



(x - Xq) 



21 



2At 



(3.4) 



approaching Eq. (3J3) for t — > 0. At longer times the bar- 
rier A/2x 1+/i comes into play preventing x from crossing 
the absorbing state x = 0. This is, however, a random 
event which can occur at an arbitrary time instant, i.e., 
the finite-time-singularity taking place at to in the deter- 
ministic case is effectively resolved in the noisy case. For 
not too large noise strength the distribution is peaked 
about the noiseless solution and vanishes for x — > 0, cor- 
responding to the absorbing state, implying the bound- 
ary condition 



P(0,t) = 



(3.5) 



In order to model a possible experimental situation the 
first-passage-time or here absorbing state distribution 
W(t) is of more direct interest (H, |29|| . First-passage 
properties in fact underlie a large class of stochastic pro- 
cesses such as diffusion limited growth, neuron dynamics, 
self-organized criticality, and stochastic resonance p8[ . 

Since P{0,t) = for all t due to the absorbing state, 
the probability that x is not reaching x = in time t 
is thus given by J Q °° P{x, t)dx, implying that the prob- 
ability —dW that x does reach x = in time t is 
— dW = — J °° dxdt(dPfdt). Consequently, the absorbing 
state distribution W(t) is determined by the expression 




W(t) = 



dP(x,t) 
Ft 



dx 



(3.6) 
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FIG. 2: In a) we show a particular realization of a;. At early 
times near xo we have random walk behavior. At longer times 
x is attracted to the absorbing state at x = 0. In b) we de- 
pict the distribution P(x, t) in the logarithmic case for /i = 0. 
For t = P(x,0) — S(x — xo), i.e., the initial condition. For 
larger t the distribution is broadened about the noiseless tra- 
jectory. P(x,t) exhibits a power law behavior for large t near 
the absorbing state. In c) we show the absorbing state dis- 
tribution W(t). For small t the distribution is exponentially 
small; for large t it displays a power law behavior for fi — 
and a stretched exponential behavior in the generic case for 
H > 0. 



In the abs ence of noise P(x, t) = S(x — x(t)) and it follows 
from Eq. ( |3.6| ) that W(t) — 5(t — to), in accordance with 
the finite time singularity at t = to. For weak noise we 
anticipate that W(t) will peak about to with vanishing 
tails for small t and large t. In Fig. || we have depicted a 
particular realization of x in the noisy case, the distribu- 
tion P(x,t) in a plot versus x and i, and the absorbing 
state distribution W(t). 

In the case of Gaussian white noise the distribution 
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P(x, t) satisfies the Fokker-Planck equation |28, B£ 



8P 
~dt 



ld_ 



—P + A— 

dx dx 



(3.7) 



in the present case subject to the boundary conditions 
P(x, 0) = S(x — .t ) and P(0, t) = 0. The absorbing state 
distribution W(t) then follows from Eq. (3.6). We note 
that the Fokker-Planck equation has the form of a conser- 
vation law dP/dt + dJ/dx = 0, defining the probability 
current J = (l/2)(dF/dx)P ~ (l/2)AdP/dx. Inserting 



Eq. (3.7) in the expression ( p.6[ ) for the distribution W(t) 
and using that J — > for x — > oo we obtain another ex- 
pression for W(t): 



W(t) = \ 



ox dx 



(3.8) 



The absorbing state distribution is thus equal to the 
probability current at the absorbing state. 



IV. WKB PHASE SPACE APPROACH 

From a stru ctural point of view the Fokker-Planck 
equation (3.7) has the form of an imaginary-time 
Schrodinger equation AdP/dt = HP, driven by the 
Hamiltonian or Liouvillian H. The noise strength A 
then plays the role of an effective Planck constant with 
P corresponding to the wavefunction. A method utiliz- 
ing a non-perturbative WKB phase space approach to a 
generic Fokker-Planck equation for extended system was 
derived in the context of the Kardar-Parizi-Zhang equa- 
tion describing interface growth |30], fnl Q . In the case 
of a single degree of freedom this method amounts to 
the eikonal approximation |^|, ^3|, |2!|, see also (3^, |34[ . 
For systems with many degrees of freedom the method 
has for example been expounded in pjjl , based on the 
functional formulation of the Langevin equation |?6|. ^| . 
In the present formulation [H], |32) the emphasis is 
on the canonical phase space analysis and the use of dy- 
namical system theory |38|, |3^| , for more details we refer 
to Appendix A. 

The weak noise WKB approximation corresponds to 
the ansatz P oc exp[— S/A]. The weight function or ac- 
tion S then to leading asymptotic order in A satisfies 
a Hamilton- Jacobi equation dS/dt + H = which in 
turn implies a princ iple of least action and Hamiltonian 
equations of motion |4C], |dj . In the present context the 
Hamiltonian has the form 



H 



1 



1 dF 
P- 



2 r 2 r dx 
yielding the equations of motion 



dx 

~dt 

dp 
~dt 



IdF 
2~dx 
1 d 2 F 
2 P lx^ 



(4.1) 

(4.2) 
(4.3) 



replacing the Langevin equation ( |l . l[ ) . The noise rj is 
then represented by th e m oment um p = dS/dx conjugate 
to x. The equations (4.2) and ( |4.3| ) determine orbits in 
a canonical phase space spanned by x and p. Since the 
system is conserved the orbits lie on the constant energy 
manifold( s) g iven b y E = H. The free energy F is given 
by Eqs. ( |2.2| ) and ( |2.3| ) and the action associated with 
an orbit from xq to x in time t has the form 



S(x -> x,t) = 



dt 



dx 



(4.4) 



According to the ansatz the probability distribution is 
then given by 



P(x, t) oc exp 



S(xq — > X, t) 



A 



(4.5) 



The zero-energy manifold E = plays an important 
role in determining the long time distributions. Insert- 
ing dF/dx = X/x 1+fJ - in Eq. (4.1) the zero-energy man- 
ifold has a submanifold structure give n by p — and 
p = \/x 1+tJ -. According to Eq. (12) the p = sub- 
manifold correspond s to the noiseless deterministic mo- 
tion given by Eq. (2.1). In Fig. || we have depicted 
the canonical phase space spanned by x and p. The 
heavy lines represent the zero-energy submanifolds p = 
and p — A/x 1+M . For E > the energy surfaces are 
equidistant for large x approaching a constant p value; 
for small x the manifold p ~ A/x 1+M for p > and 
p — -AEx 1+tl /X for p < 0. For E < the energy 
surfaces are confined between the zero-energy subman- 
ifolds; the manifolds approach (x,p) = (0,0) according 
to p ~ 4\E\x 1+,1 /\ and for large p as p ~ \/x 1+tl . 
For E — > — oo the orbits approach the positive p half- 
axis. The arrows indicate the direction of motion on 
the manifolds. The dashed line indicates a nullcline 
(dx/dt = 0) passing through the hyperbolic fixed point 
(x,p) = (oo, 0). In the long time limit the orbit from Xq 
to x converges towards the zero-energy submanifolds. 



V. RANDOM WALK AND LONG TIME TAILS 

The weak noise phase space approach reviewed above 
affords a simple derivation of the asymptotic long 
time behavior of the distributions for the finite-time- 
singularity problem. In order to derive the transition 



probability P(x, t) according to Eq. (4.5) we simply have 
to identify the relevant orbit in phase space from xq to 
x which at long times passes close to the zero energy 
manifolds. 



A. The random walk case 

It is instructive first to consider the case A = 0. Here 
the finite time singularity at x = is absent, there is no 
absorbing state and the Langevin equation ( |1.1| ) takes 
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corresponding to constant values of p. The bold line for p — 
is the zero-energy manifold. 



FIG. 3: We show the topology of phase space. The bold lines 
indicate the zero energy submanifolds. We show representa- 
tive orbits for E > and E < 0. The dashed line indicates a 
nullcline (dx/dt = 0) to the saddle point (x,p) = (oo,0). At 
long times the orbit from xo to x migrates to the zero energy 
submanifolds. 



This behavior is in accordance with the general discus- 
sion in Section III and is graphically depicted in Fig. ||. 
The phase space topology in the random walk case is 
shown in Fig. 0. 



the form dx/dt — rj(t), describing random walk on the 
whole axis. The Hamiltonia n is given by H = (l/2)p 2 
and the equations of motion (4.2) and (4.3) have the form 
dx/dt — p, dp/dt — with solution s p = po and x = xq + 
Pot. Inserting dx/dt and H in Eq. (4.4) for the action we 
obtainS' = (1/2) J^dtp 2 = {l/2)p 2 l t and finally using Eq. 



(4.5) the Gaussian distribution ( p.4[ ) for random walk. 

We note that in order to obtain the correct limit of 
the finite-time-singularity problem we must incorporate 
the absorbing state condition at x = 0. This is achieved 
by using the method of mirrors p3[ and considering for 
P(x, t) the linear combination, 



P(x,t) = {2nAt) 
(x - x ) 2 



-1/2 



exp 



2At 



exp 



(x + Xo) 



2At 



, (5.4) 



in the half space x > 0. This distribution is a solution of 
the Fokker Planck equation and vanishes for x = 0. For 
small x it behaves linearly with x, 

P(x,t) = (2/7r) 1/2 x(At)- 3/2 iEoexp(-.^/2Af) . (5.2) 



Using Eq. ( p.8| ) we readily obtain the well-known random 
walk result 

W{t) = (2/Tr) 1 ^x Q (At)-^ 2 exp(-x 2 /2Ai) . (5.3) 

For small t the distribution vanishes exponentially. It 
displays a maximum at to = x^/3A and falls off alge- 
braically as t~ a for large t with scaling exponent a = 3/2. 



B. The absorbing state case 



For zero noise and for A = we have from Eq. (4.1) 
x = xo at all times, whereas the solution in the case of a 
finite-time-singularity is attracted to the absorbing state 
at time to- In the presence of noise the attraction gives 
rise to a change of the form of the absorbing state distri- 
bution W(t) from a delta function peak to a broadened 
peak. For large A the distribution shows a maximum 
about to; for intermediate values of A the maximum is 
between to and the random walk value Xq/3A. Because 
of the attraction to the absorbing state we also obtain a 
faster long time fall-off and thus a positive correction to 
the random walk exponent a = 3/2, depending on the 
strength A. 

The plot in Fig. [| permits a simple qualitative dis- 
cussion of the finite-time-singularity phase space phe- 
nomenology. Firstly, a short time orbit from xq to x 
corresponds to large value of \p\. In this region and for 
not too small x the phase space topology in Fig. || is 
similar to the random walk case depicted in Fig. ^ and 
we infer u nbia sed ra ndo m walk behavior, yielding the ex- 
pression (5.1) and (|5.3|). Secondly, searching for longer 
time orbits from £o to x we must choose smaller p and we 
move into the region of phase space for negative energy 
where the finite-time-singularity dominates the topology, 
as shown in Fig. ||[ In the limit of long times the orbit 
approaches asymptotically the zero-energy submanifolds. 
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1. The logarithmic case /i = 

In the logarithm ic c ase /i = the zero energy con- 
dition using Eq. (11) and Eq. (2.2) yields the rela- 
tionship p = X/x, corresponding to the hyperbolic man- 
ifold; note that the p — manifold corresponds to de- 
terministic motion and yields S = 0. Setting p = X/x 
and H = in the action in Eq. (4.4) we then have 
S = J dt pdx/dt — J x dx X/x = A l og a: . Moreover, 
for p — X/x the equation of motion (4.2) reduces to 
dx/dt = X/2x with the growing solution x 2 = Xt. In 
the long time limit where the orbit is close to the zero- 
energy manifold we thus obtain S ~ (A/2)logt, yield- 
ing according to Eq. ( |4.5| ) the power law distribution 

P(x,t) oc t~^s. Owing to the absorbing state the dis- 
tribution P(x, t) must vanish for x — > 0. As discussed 
in Appendix ^ this limit is reproduced by pushing the 
WKB approximation in Section IV to next order in A, 
yielding the correction — A log a; to S. Finally, we obtain 
in the weak noise- long time limit the distribution 



The expressions (5.6) and ( |5.7| ) show that in the case of a 
generic finite-time singularity characterized by the index 
ll the power law behavior of P(x, t) and W(t) is altered to 
a stretched exponential behavior depending on ll. In the 
limit ll — ► we obtain the power law behavior. We note 
again that the WKB approximation is unable to produce 
the random walk prefactor i~ 3 / 2 for A = and the peak 
of W(t) about t . 



VI. SOLUTION OF THE FOKKER-PLANCK 
EQUATION 

In this section we return to the Fokker-Planck equa- 
tion ([O]) and present exact expressions for the transition 
probability P(x, t) and the absorbing state distribution 
W(t) in the the logarithmic case ll = 0. We have summa- 
rized key points in the derivation here and defer details 
to Appendix IbI. 



P(x, t) oc xt 



(5.4) 



Moreover, applying Eq. ( |3.8| ) we deduce the weak noise- 
long time absorbing state distribution 



W(t) oc r 



(5.5) 



with sca ling exponent a = A/2A The expressions (5.4) 
and (5.5) show that the finite-time-singularity or, equiv- 
alently, absorbing state attracts the random walker and 
increase the fall-off exponent a. We note that the WKB 
approximation to leading order in A fails to retrieve the 
random walk exponent 3/2 in the limit A — > and the 
maximum of W(t) about to. 



A. Quantum particle in a repulsive potential 

In the logarithmic case the Fokker-Planck equation as- 
sumes the form 



dP _ A d 2 P X OP A 
"dt ~ ~2~dx 2 + 2x ~dx ~ 2^ 



P 



(6.1) 



Removing the first order term by means of the gauge 
transformation exp(/i) = x~ x / 2A we have 



d 

- A— [exp(-/i)P] = H [exp(-/i)P] 
where the Hamiltonian H is given by 



(6.2) 



2. The generic case /j, > 

In the generic case li > the zero-energy mani- 
fold implies the constraint p — A/x 1+M and we ob- 
tain similar to the logarithmic case above the action 
S = fdt pdx/dt = J x dx X/x 1+tl = -(\/(jl)x-i* + const, 
and the equation of motion, dx/dt = A/2x 1+M , on the 
zero-energy manifold with solution x M+2 = (1 + (i/2)Xt. 
In the long time limit we thus find the action S = 
-(A//i)((l + f i/2)Xt)-»/( 2+ ^ + X/ll. The second order 
correction to S, evaluated in Appendix is given by 
— A(l + //) logcc and we obtain the weak noise-long time 
distribution 



P(x,t) oc 
exp 



A 

All 



(( 



1 



^Xt) 2+ "-i 



and absorbing state distribution 

" A 



W{t) - exp 



All 



(5.6) 



(5.7) 



H 



1 A2 9 2 

-A 2 

2 dx 2 



X 2 r 

8 



1 



2A 



(6.3) 



This Hamiltonian describes the motion of a unit mass 
quantum particle in one dimension subject to a centrifu- 
gal barrier of strength (A 2 /8)(l + 2A/A). For A = the 
barrier is absent and the particle can move over the whole 
axis; this case corresponds to ordinary random walk [p3| . 
For A ^ the particle cannot cross the barrier and is con- 
fined to either half space; this corresponds to the case of a 
finite-time-singularity subject to noise and an absorbing 
state at x = 0. 

The Fokker-Planck equation ( |6.2[ ) has the form of an 
imaginary time Schrodinger equation with Planck con- 
stant A for the wavefunction exp(— h)P and is readily 
analyzed in terms of Bessel functions |42]| . Incorporat- 
ing the initial condition P(x, 0) = 5(x — Xo) by defining 
P(x, t) — > P(x, t)8(t) we obtain the inhomogeneous dif- 
ferential equation 



dP(x, t) 

Ft 



5{x — xo)S(t) 



? h He- h P(x,t) , (6.4) 
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for the determination of P(x,t). On the positive x and 
k axis the wavefunctions ipk( x ) — (kx) 1 / 2 J\/2+\/2A{kx) 
form, according to the Fourier-Bessel transform 0], an 
orthonormal and complete set satisfying the eigenvalue 
equation Ht/jk {%) — (A 2 /2)k 2 ijjk(x). Expanding the right 
hand side of Eq. ( |6.4[) on the set ipk and using a well- 
known identity for Bessel functions Q, [|5| we obtain for 
the probability distribution P(x, t) the following closed 
expression: 



C. The distribution W(t) 



The absorbing state distribution W{t) in Eq. (6.6) 
vanishes exponentially for small t. For large t the dis- 
tribution shows a power law behavior with exponent 
a = 3/2 + A/2A. In the deterministic limit A = we 
have W(t) = S(t — to), where to = Xq/X. For small A the 
exponent approaches A/2A in accordance with Eq. ( |5.5| ). 
The distribution has a maximum at 



A ,1 (^ + *0> 

p( xt ) = — r n, x {— -) . 6.5) 

„2s-2 At 2 + \ At J y ' 



tu 



X, 



(I 



Here l v is the Bessel function of imaginary argument, 

i u ( z ) = {-iyjSz)Q>- 

By means of Eq. ( |3.8| ) we moreover deduce the absorb- 
ing state distribution 



W(t) 



2.7.:, 



1 + A/A 




T(l/2 + A/2A) 



(6.6) 



Similar expressions have also been derived in the context 
of the XY model 151. 



B. The distribution P(x,t) 



The expression (5.5) provides the complete solution of 
the finite-time-singularity problem for \i = 0. The ex- 
pression is discussed in more detail in Appendix [^. For 
t = we have P(x, 0) = S(x — x ) in accordance with the 
initial condition ( |3.3[ ) . For small t we obtain the random 
walk result P(xJ)= exp(-(x - x ) 2 /2At)/(2nAt) 1 / 2 
in accordance with Eq. (3.4). For A = we have 
P(x,t) = (exp(-(x-x ) 2 /2At) +x a -> -x a )/(2nAt) 1 ^ 2 
in agreement with Eq. (5.1) for random walk with an 
absorbing wall at x = 0. 

For long times and x close to the absorbing state x = 
we obtain the asymptotic form 



P(x,t)<x—± 



2A) 



-(2Ai)- 



(6.7) 



For small x the distribution vanishes linearly due to the 
absorbing state. For large t the distribution exhibits a 
power law behavior with scaling exponent a = 3/2 + 
A/2A. In the weak noise limi t th e distribution is peaked 
about the noiseless solution (2.4) . Fo r A = we obtain 
the random walk result in Eq. ( |5.2D and we note that 
the finite-time-singularity or absorbing state lead to an 
increase of the scaling exponent and thus a faster fall-off 
in time. In the weak noise limit the scaling exponent 
approaches A/ 2 A in agreement with the WKB analysis 
m Section 0. In Fig. [2j we have depicted the distribution 
P(x,t). 



3A + A 



(6.8) 



For A = we have t max — > to and for A = the random 
walk result i max = xq/3A. For large coupling strength 
tmax 0. Expanding W(t) about i max we obtain the 
Gaussian distribution 



W{t) cx e -(*-W ) 2 /- 2 i 
characterized by the mean square width 

2 _ iAx 

° ~ (3A + A) 3 ' 



(6.9) 



(6.10) 



Since W(t) falls off as a power of t only a finite number 
of moments (t n ) = Jt n W(t) dt exists. For (2n-l)A < A 
we have 



(n = n 



(6.11) 



XXVA-(2p-l)A 
The distribution W{t) is shown in Fig. |2[ 

VII. SUMMARY AND CONCLUSION 

In this paper we have addressed the problem of the in- 
fluence of white Gaussian noise of strength A on a generic 
finite-time-singularity of strength A, characterized by the 
exponent [i. We have for simplicity considered only a sin- 
gle degree of freedom. We have found that in the case 
of a logarithmic sink in the free energy driving the vari- 
able, corresponding to a square root singularity, the first- 
passage- time or absorbing state distribution W(t) dis- 
plays a peak about the finite-time-singularity and a long 
time power law tail oc i~ Q , characterized by the scaling 
exponent a = 3/2 + A/2A. The exponent is nonuniversal 
and depends on the ratio between the singularity strength 
A and the noise strength A. In the case where the noise 
originates from a thermal environment at temperature T 
we have A cx T and the scaling exponent depends on the 
temperature, a = 3/2 + const. /T. 

In the generic case of a finite-time-singularity charac- 
terized by the exponent fi > the weak noise WKB ap- 
proach shows that the power law tail for fi = is changed 
to a stretched exponential with a slower fall-off. 

To the extent that the character of a finite-time- 
singularity in the vicinity of threshold can be modeled 
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with a single degree of freedom the present study should 
hold as regard the influence of noise on the time distribu- 
tion. We note in particular that in the case of a thermal 
environment at temperature T the change of the scaling 
exponent becomes large in the limit of low temperatures 
as the distribution narrows around the noiseless thresh- 
old time. 

The present study also suggests generalizations to the 
case of damping and to the case of several coupled vari- 
able subject to a finite-time-singularity. 



Acknowledgments 

Discussions with D. Diakunov, A. Brandenburg, J. 
Krug, T. Bohr, M. H. Jensen, A. Svane and K. M0lmer 
are gratefully acknowledged. We also thank A. Bray, D. 
Sornette, and R. Mantegna for drawing our attention to 
similar results within the context of the critical dynamics 
of the XY model, modeling with finite-time-singularities 
in econophysics and related fields, and long-range corre- 
lated stationary Markov processes, respectively. 



APPENDIX A: THE PHASE SPACE METHOD 

The weak noise WKB approximation applied to the 
Fokker-Planck equation is well-documented |2^, |2^, 
|34| |. Here we review this method with emphasis on a 
canonical phase space approach which we have found use- 
ful in discussing the pattern formation and scaling in the 
noisy Burgers equation [po|, p[ [to], p| p|. We also note 
that the approach follows from a saddle point approxi- 
mation to the functional Martin-Siggia-Rose approach to 
nonlinear Langevin equations fl36L t37L f49, M . 



1. To leading order A 

Taking as our starting point a generic Langevin equa- 
tion for one degree of freedom x driven by Gaussian white 
noise, 

^- = -l G {x)+r 1 {t) , ( VV )(t)=5(t), (Al) 

the associated Fokker-Planck equation for the distribu- 
tion P(x, t) takes the form 



dP _ ld_ 
~dt ~ 2dx 



GP + A 



dP 

dx 



Applying like in quantum mechanics 
proximation 



P{x, t) (x exp 



S(x,t) 



(A2) 

the WKB ap- 
(A3) 



and expanding the action S in powers of the noise 

strength A, S = Sq + AS±, Sq satisfies the Hamilton- 
Jacobi equation jlO| 

f +*<*«>-<>. » = §• (A4) 

implying a principle of least action pcjj . The action So 
thus has the form 

S (x,t) = J dt (p^-h) , (A5) 



with Hamiltonian given by 

H=\ P {p-G) , 

implying the Hamilton equations of motion 
dx 1 



dt 
dp 



2 

I ^ 

2 P dx 



G+p , 



(A6) 

(A7) 
(A8) 



The deterministic coupled equations (A7) an d (A£) re- 
place the stochastic Langevin equation (1.1) for weak 
noise. The noise rj is replaced by the canonic ally con- 
jugate momentum p, which by means of Eqs. ( |A3| ) and 
( A4 ) is given in terms of the distribution P 



p = —Ad\ogP/dx 



(A9) 



The equations of motion define orbits in a (x,p) phase 
space lying on the constant energy surfaces E — H and 
the general discussion of the original stochastic problem 
is replaced by an analysis of the phase space topology. 
The prescription for deriving the distribution to leading 
order in A thus amounts to i) solve the equations of mo- 
tion ( A 7 ) and ( A8 ) for an orbit from an initial value Xq 
to a final value x reached in the time span t, p being a 
slaved variable, ii) evaluate the action S associated with 
an orbit according to Eq. (A5), and, finally, hi) derive 
the tran siti on probability from xq to x in time t using the 
ansatz ( A3 ) . The zero-energy manifold here play an im- 
portant role in determining the long time distributions. 
In the limit t — > oo a given finite time orbit from xq to x 
thus converges to the zero-energy manifold. 

In the case of an overdamped harmonic oscillator de- 
scribed by the Langevin equation dx/dt = —ujx + r] the 
phase space analysis was carried out in |3(J . In the case 
of random walk given by dx/dt = rj the analysis is per- 
formed in Section V A . 



a. Finite-time- singularity case 

In the generic finite-time-singularity case, G = 
A/x 1+M , and we have the Hamiltonian H — (l/2)p(p — 
A/x 1+M ), yielding the equations of motion dx/dt = 
-A/2x 1+ ^ +p and dp/dt = -A(l + fi)p/2x 2+t ". These 
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equations are, however, not particularly tractable and 
we therefore only consider the logarithmic case for ji = 0. 
Here G — X/x and we obtain the Hamiltonian 



H = \v 1 /' 
and the equations of motion 



dx 


X 


~dt ~ 




dp 


A p 


~dt ~ 


2x 2 ' 



(A10) 

(All) 
(A12) 



These equations have a hyperbolic fixed point at 
(x,p) = (oo,0). The nullcline dx/dt — to the sad- 
dle point is given by p — X/2x, indicated in Fig. || by 
the dashed line. The zero-energy manifolds are given by 
p = and p = X/x. The conserved energy H = E pro- 
vides the first constant of integration. Solving for x we 
have x = X/{p 2 - 2E). For E > p ->■ ±{2E) 1 / 2 for 
x — > oo as indicated in Fig. |^, yielding in that limit the 
random walk phase space topology depicted in Fig. ||. 
For E < x = X/(p 2 + 2\E\) exhibiting a maximum at 
the nullcline in a plot of x versus p as shown for two rep- 
resentative orbits in Fig. |3[ Using energy conservation to 
solve the equation of motion for p and subsequently for x 
we obtain the solutions x 2 = (t + ii)(A + 2E(t + ti)) and 
p 2 = 2E + X/(t + ii), where t\ is the second constant of 
integration. A specific orbit from xq to x\ in time t thus 
determines the constants E and t\ ; the momentum p be- 
comes a slaved variable, and the action evaluated along 
the orbit yields the distribution. 

Considering as final value the absorbing state x\ = 0, 
the long time orbits lie in the negative energy region and 
eliminating E we obtain the solution, < t' < t 



= (xl (1 - t'/t) + At') (1 - t'/t) 



(A13) 



The energy is given by \E\ = (Xt — xfy/t 2 and we note 
that the energy approaches zero in the long time limit, 
i.e., the orbit from x to x — migrates to the zero- 
energy manifold. Finally, for the action associated with 
the orbit we obtain 



S = 



1 



Alog — - 1 



yielding the long time distribution 
P{x ^0,t)(x 



(A14) 



(A15) 



in accordance with the expression (5.4). 

Alternatively, elim inatin g the momentum p the equa- 
tions of motion ( All ) and ( A12) reduce to a second order 
equation for x, d l xj ' dt 2 — —dV/dx, describing the mo- 
tion of a particle of unit mass in the attractive potential 
V{x) = —(l/8)X 2 /x 2 . It then follows by simple quadra- 
ture that all direct orbits to the absorbing state x = 
take a finite time, whereas the traversal time of negative- 
energy orbits with a turning point diverges in the limit 
\E\ — > 0; this is in accordance with the phase space be- 
havior shown in Fig ^. 



2. Next leading order in A 

The next leading order in A is obtained from Si which 
by insertion satisfies the equation of motion 



dSi 
~~dt 



dSo 
dx 



G\ dSi 
~2)~dx~ 



1 d 



G 



dS Q 



(A16) 



2 da; V dx 
where dSo/dx is obtained from the first order solution. 

a. Random walk case 

From the random walk case discussed in Sectio n | VA| 
we have So = ( x ~ %o ) 2 /2t. Consequently, Eq. (|A16Q 
takes the form 



dSi 
~~dt 



X — Xq \ dSi 1 



t 



t I dx 2 
with a particular time-dependent solution 

Si = ~log|t| , 



(A17) 



(A18) 



yielding S = (x — x ) 2 /t + (A/2) log \t\ and the Gaussian 
distribution 



P(x,t) oc \t\- x /' 2 



exp 



(x - x ) 2 
At 



(A19) 



As in the quantum case p2f the next leading correction 
yields the normalization factor jij" 1 / 2 . 



b. Finite-time- singularity case 

In the finite-time-singularity case G = A/x 1+M and 
from above dSo/dx — 0. We then obtain inserting in 
Eq. (|A16|) 



dSi 

~dT 



X dSi 1 + n A 



2x 1+ ^ dx 2 x 2 +v 
with a particular space-dependent solution 

Si = -{l + n)\ogx , 

giving rise to the factor x 1+M in Eq. (|5~ 



(A20) 



(A21) 



APPENDIX B: THE FOKKER PLANCK 
EQUATION 



Here we discuss the Fokker-Planck equation (3.1) in 
the logarithmic case /i = in more detail. Applying the 
gauge transformation, exp(/i) — x ~ A / 2A , and incorporat- 
ing the boundary condition ( |3.3| ), P(x, 0) = 6(x — xq), 
we obtain the inhomogeneous differential equation ( |6.4| ), 
with Hamiltonian H given by Eq. ( ]6.3D corresponding to 
the motion of a quantum particle subject to a centrifugal 
barrier cx 1/x 2 . 
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Exact solution 



2. Random walk, short time, and long time limits 



The right hand side of Eq. (|6.l | ) has the same form 
as the standard Bessel equation |ff4| ft5| . Noting also 
the analogy to the quantum case of particle motion in 
spherical coordinates [42] it follows that 



ipk(x) = {kx) 1/2 Z y {kx) 



(Bl) 



where Z v (kx) is a solution of the Bessel equation, satisfies 
the eigenvalue equation 



Hip k (x) = k 2 tp k (x) , 



(B2) 



for v = ±(1/2 + A/2 A). The Bessel function of the first 
kind J v (kx) satisfies the absorbing state boundary condi- 
tion J v — > for x — + and the completeness and orthog- 
onality of ipkjx) follow from the Fourier-Bessel integral 
representation Hj 

/>00 />OG 

f(r) = dk dr'kJ v {kr)J v {kr')f{r') , (B3) 
Jo Jo 

valid for v > 1/2. We proceed by Fourier transforming 



P(x,t) 



2^ 



(B4) 



and subsequently expanding p w (x)x A ^ 2A and S(x — Xq) 
on the eigenfunctions ipk{x), 

/>OC 

p w (x)x A/2A = / dkMx)P^k , (B5) 



S(x-x ) = dkipk(x)i/Jk(xo) , (B6) 
Jo 



yielding the expansion coefficients 

^o /2A ^fc(a:o) 
-iw + Afc 2 /2 ' 



(B7) 



Finally, integrating over lo in the lower half plane and 
picking up contributions from the branch cut we obtain 

P(x,t) = 

dk fce- Afc2 */ 2 (xx ) 1/2 (z/zo) A/2A x 

Jl/2+X/2A(kx)J 1/2+ X/2A(kxo) ■ (B8) 

This integral can be reduced further using the identity 
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dx xe p x J u (ax)J v {(ix) 



l_ e io?^)l^ Iv{a(}l2(?) 



(B9) 



yielding Eq. (p.q) in Section VI 



In the random walk case for A = we obtai n us - 
ing Ii/ 2 — 2(l/2irx) 1 / 2 sinhcc, the expression (|5.l|). 

In the short time limit t <C xxq/A, using I v {x) oc 

(1/27TX) 1 / 2 exp(x) for x — > oo [||| we obtain Eq. ( |3.4| ) 
and for i = the boundary condition (3.3). 

In the long time limit t xxq/A usin g Iy (x) oc 
(x/2y/r(v + l) for x -> @ we obtain Eq. Q Using 
Eq. (3.S) and T(z + 1) = zTJz) we finally obtain the 
absorbing state distribution (p.6|) . 

The moments of 1U(£) are easily worked out. Using 
T(z) = J °° t z - x e- 1 dt g H we have 

(* n ) = / * - (ff )" r( | ( tf r } n) '( bi °) 

or fu rther reduced for (2n — 1)A < A the expression 



Weak noise limit 



In the limit A — > the dis tribu tion P(x,t) is centered 
about th e no iseless solution (2.4). In terms of the exact 
solution (6.5) this is a singular limit since both order and 
argument in I v {x) diverge. Using the spectral represen- 
tation @, |4| 



I v (z) = 



{z/2Y 



2Y C 

1 x I cosh(x cos 9) sin 2l/ dO , 
2> T (2' Jo 

(Bll) 



introducing the variable u according to 

Af 



sinh u 



2xxo 



defining 



/±(( ,, = 1 „ gsta , + I( 1± i^- 



and using 



r i 



A 



, T I - ) « Trx/2e-^ , 
2A/ \2 V 2A 



i \ 2+2A 



(B12) 



(B13) 



,(B14) 



for small A we obtain by insertion in Eq. (6.5) 
P(x,t) 



1 1 \ /„2\ 2 + 23 



4irV2 xo A V At 



u \ „-(i 2 +i5)/2Ai 



d(9 — — 

i , sinh u 



(B15) 



The expression (B15) for P(x,t) is directly amenable to 



an asymptotic analysis for A — > by means of Laplace's 
method IplJ. For small A the main contributions to the 
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integral originate from the maxima of f+(0) and /-(#)■ 
The two maxima in the interval < 8 < 7r are given 
by cos6>± = ±exp(— u), yielding f±(9±) = — cothw and 
f±(0±) = (l/2)(log(l - e- 2u ) + cothu). Performing the 
Gaussian integrals about the maxima we thus obtain the 
asymptotic result valid for small A and fixed u, i.e., fixed 
x/t. 



P(x,t) 
(1-e- 



A 



2vrA 



(sinh2u)2 



exp 



x \XtJ 
2At 



1 , A 
2+ TK 



(B16) 



where F(x, t) = x 2 + x 2 , — Xt cothw or by insertion 

F(x,t) = x 2 +xl - [(2xx ) 2 + (Xt) 2 ] 1/2 . (B17) 



In the short time limit At <C 2xxq u is small and 
sinhw ~ u, i.e., u ~ Xt/2xxo and we obtain P{x,t) ~ 
(l/27rAt) 1 / 2 exp(-(x - x ) 2 /2Ai), yielding 5{x - x ) 
for t = 0. The weak coupling limit A — > for 
fixed x and t is also consistent; we obtain P(x, i) ~ 
(l/27rAt) 1 / 2 exp(-(x - x ) 2 /2Ai). For weak noise the 
peak of the distribution is determined by the condition 
F = 0, yielding x = (x 2 , — At) 1 / 2 and the peak thu s follows 
the noiseless finite-time-singularity solution (2.4). 
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